Loading...
 

Algorytm eliminacji Gaussa z pivotingiem

Zacznijmy ponownie eliminacje Gaussa dla naszego przykładu, tym razem z zaminą wierszy (zwaną również "pivotingiem" w dosłownym tłumaczeniu z języka angielskiego).


\( \begin{bmatrix} {\color{blue}0.05} & 0.10833 & 0.00833 & 0 & 0 \\ {\color{red}0.10833} & 0.5 & 0.21666 & 0.00833 & 0 \\ {\color{blue}0.00833} & 0.21666 & 0.55 & 0.21666 & 0.00833 \\ {\color{blue}0} & 0.00833 & 0.21666 & 0.5 & 0.10833 \\ {\color{blue}0} & 0 & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix}1 \\ 1 \\ 1 \\ 1 \\ 1 \end{bmatrix} \)

Rozpoczynamy od operacji zamiany wierszy. W pierwszej kolumnie (zaznaczonej na niebiesko i czerwono) szukamy największej wartości na moduł. Jest nią wyraz \( A_{2,1}=0.10833 \) (zaznaczony na czerwono). Zamieniamy więc wiersz pierwszy z drugim.

\( \begin{bmatrix} {\color{red}0.10833} & 0.5 & 0.21666 & 0.00833 & 0 \\ {\color{blue}0.05} & 0.10833 & 0.00833 & 0 & 0 \\ {\color{blue}0.00833} & 0.21666 & 0.55 & 0.21666 & 0.00833 \\ {\color{blue}0} & 0.00833 & 0.21666 & 0.5 & 0.10833 \\ {\color{blue}0} & 0 & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix}1 \\ 1 \\ 1 \\ 1 \\ 1 \end{bmatrix} \)

Teraz, wykonujemy krok "normalnego" algorytmu eliminacji Gaussa.
Skalujemy pierwszy wiersz tak żeby dostać jedynkę na przekątnej
\( 1^{st}=1^{st}A_{1,1}=1^{st}\frac{1}{0.10833}=1^{st}9.231053 \)
\( \begin{bmatrix} 0.10833{\color{red}*9.23105} & 0.5{\color{red}*9.23105} & 0.21666{\color{red}*9.23105} & 0.00833{\color{red}*9.23105} & 0{\color{red}*9.23105} \\ 0.05 & 0.10833 & 0.00833 & 0 & 0 \\ 0.00833 & 0.21666 & 0.55 & 0.21666 & 0.00833 \\ 0 & 0.00833 & 0.21666 & 0.5 & 0.10833 \\ 0 & 0 & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 1{\color{red}*9.23105} \\ 1 \\ 1 \\ 1 \\ 1 \end{bmatrix} \)
\( \begin{bmatrix} {\color{red}1} & {\color{red}4,61553} & {\color{red}2.00006} & {\color{red}0,07693} & {\color{red}0} \\ 0.05 & 0.10833 & 0.00833 & 0 & 0 \\ 0.00833 & 0.21666 & 0.55 & 0.21666 & 0.00833 \\ 0 & 0.00833 & 0.21666 & 0.5 & 0.10833 \\ 0 & 0 & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix}{\color{red}9.2311} \\ 1 \\ 1 \\ 1 \\ 1 \end{bmatrix} \)

Następnie odejmujemy pierwszy wiersz od drugiego tak by dostać zero w pozycji \( A_{2,1} \)
\( 2^{nd}=2^{nd}-1^{st}A_{2,1}=2^{nd}-1^{st}0.05 \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ {\color{red}0} & 0.10833{\color{red}-4.61553*0.05} & 0.00833{\color{red}-2*0.05} & 0{\color{red}-0.07693*0.05} & 0 \\ 0.00833 & 0.21666 & 0.55 & 0.21666 & 0.00833 \\ 0 & 0.00833 & 0.21666 & 0.5 & 0.10833 \\ 0 & 0 & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \\ \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix}9.23105 \\ 1{\color{red}-9.23105*0.05} \\ 1 \\ 1 \\ 1 \end{bmatrix} \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ {\color{red}0} & {\color{red}-0.12244} & {\color{red}-0.09167} & {\color{red}-0,00385} & 0 \\ 0.00833 & 0.21666 & 0.55 & 0.21666 & 0.00833 \\ 0 & 0.00833 & 0.21666 & 0.5 & 0.10833 \\ 0 & 0 & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ {\color{red}0.53844} \\ 1 \\ 1 \\ 1 \end{bmatrix} \)
Odejmujemy teraz pierwszy wiersz od trzeciego tak by dostać zero w pozycji \( A_{3,1} \)
\( 3^{rd}=3^{rd}-1^{st}A_{3,1}=3^{rd}-1^{st}0.00833 \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & -0.12244 & -0.09167 & -0,00385 & 0 \\ {\color{red}0} & 0.21666{\color{red}-4.61552*0.00833} & 0.55{\color{red}-2*0.00833} & 0.21666 {\color{red}-0.07689*0.00833} & 0.00833{\color{red}-0} \\ 0 & 0.00833 & 0.21666 & 0.5 & 0.10833 \\ 0 & 0 & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \\ \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 0.53844 \\ 1{\color{red}-9.23105*0.00833} \\ 1 \\ 1 \end{bmatrix} \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & -0.12244 & -0.09167 & -0,00385 & 0 \\ {\color{red}0} & {\color{red}0.17821} & {\color{red}0.53334} & {\color{red}0.21603} & {\color{red}0.00833} \\ 0 & 0.00833 & 0.21666 & 0.5 & 0.10833 \\ 0 & 0 & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 0.53844 \\ {\color{red}0.92310} \\ 1 \\ 1 \end{bmatrix} \)
Ponownie wykonujemy teraz operacje zamiany wierszy.
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & {\color{blue}-0.12244} & -0.09167 & -0,00385 & 0 \\ 0 & {\color{red}0.17821} & 0.53334 & 0.21603 & 0.00833 \\ 0 & {\color{blue}0.00833} & 0.21666 & 0.5 & 0.10833 \\ 0 & {\color{blue}0} & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix}9.23105 \\ 0.53844 \\ 0.92310 \\ 1 \\ 1 \end{bmatrix} \)
W drugiej kolumnie od przekątnej w dół (zaznaczonej na niebiesko i czerwono) szukamy największej wartości na moduł. Jest nią wyraz $A(3,1)=0.17821$ (zaznaczony na czerwono) ponieważ \( |0.17821|>|-0.12244|=0.12244 \). Zamieniamy więc wiersz drugi z trzecim.
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & {\color{red}0.17821} & 0.53334 & 0.21603 & 0.00833 \\ 0 & {\color{blue}-0.12244} & -0.09167 & -0,00385 & 0 \\ 0 & {\color{blue}0.00833} & 0.21666 & 0.5 & 0.10833 \\ 0 & {\color{blue}0} & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 0.92310 \\ 0.53844 \\ 1 \\ 1 \end{bmatrix} \)
Teraz, skalujemy nowy drugi wiersz tak by na przekątnej dostać jedynkę
\( 2^{nd}=2^{nd}\frac{1}{A_{2,2}}=2^{nd}\frac{1}{0.17821}=2^{nd}5.61135 \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & 0.17821{\color{red}*5.61135} & 0.53334{\color{red}*5.61135} & 0.21603{\color{red}*5.61135} & 0.00833{\color{red}*5.61135} \\ 0 & -0.12244 & -0.09167 & -0,00385 & 0 \\ 0 & 0.00833 & 0.21666 & 0.5 & 0.10833 \\ 0 & 0 & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \\ \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 0.92310{\color{red}*5.61135} \\ 0.53844 \\ 1 \\ 1 \end{bmatrix} \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & {\color{red}1} & {\color{red}2.99275} & {\color{red}1.21222} & {\color{red}0.046742} \\ 0 & -0.12244 & -0.09167 & -0,00385 & 0 \\ 0 & 0.00833 & 0.21666 & 0.5 & 0.10833 \\ 0 & 0 & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ {\color{red}5.17983} \\ 0.53844 \\ 1 \\ 1 \end{bmatrix} \)
Odejmujemy teraz drugi wiersz od trzeciego przeskalowny tak by dostać zero w pozycji \( A_{3,2}=-0.12244 \)
\( 3^{rd}=3^{rd}-2^{nd}A_{3,2}=3^{rd}-2^{nd}(-0.12244) \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & 1 & 2.99287 & 1.21226 & 0.04676 \\ 0 & {\color{red}0} & -0.09167{\color{red}+2.99275*0.12244} & -0,00385{\color{red}+1.21222*0.12244} & 0{\color{red}+0.046742*0.12244} \\ 0 & 0.00833 & 0.21666 & 0.5 & 0.10833 \\ 0 & 0 & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \\ \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 5.17983 \\ 0.53844{\color{red}+5.17983*0.12244} \\ 1 \\ 1 \end{bmatrix} \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & 1 & 2.99287 & 1.21226 & 0.04676 \\ 0 & {\color{red}0} & {\color{red}0.27479} & {\color{red}0.14459} & {\color{red}0.00573} \\ 0 & 0.00833 & 0.21666 & 0.5 & 0.10833 \\ 0 & 0 & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 5.1800 \\ {\color{red}1.1727} \\ 1 \\ 1 \end{bmatrix} \)
Podobnie odejmujemy drugi wiersz od czwartego
\( 4^{th}=4^{th}-2^{nd}A_{(4,2}=4^{th}-2^{nd}0.00833 \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & 1 & 2.99287 & 1.21226 & 0.04676 \\ 0 & 0 & 0.27476 & 0.14457 & 0.005723 \\ 0 & {\color{red}0} & 0.21666{\color{red}-2.99275*0.00833} & 0.5{\color{red}-1.21210*0.00833} & 0.10833{\color{red}-0.046742*0.00833} \\ 0 & 0 & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \\ \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 5.17983 \\ 1.17265 \\ 1{\color{red}-5.17983*0.00833} \\ 1 \end{bmatrix} \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & 1 & 2.99287 & 1.21226 & 0.04676 \\ 0 & 0 & 0.27476 & 0.14457 & 0.005723 \\ 0 & {\color{red}0} & {\color{red}0.19173} & {\color{red}0.48990} & {\color{red}0.10794} \\ 0 & 0 & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix}9.23105 \\ 5.17999 \\ 1.17270 \\ {\color{red}0.95684} \\ 1 \end{bmatrix} \)
Została nam podmacierz 3 na 3 z prawego dolnego rogu, dla której kontynuujemy eliminację Gaussa z zamianą wierszy. Szukamy najpierw maksymalnej (na moduł) wartości w trzeciej kolumnie od przekątnej w dół (zaznaczonej na czerwono i niebiesko). Jest nią wartość z przekątnej \( A_{3,3} \) zaznaczona na czerwono, nie musimy więc zamieniać wierszy.
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & 1 & 2.99287 & 1.21226 & 0.04676 \\ 0 & 0 & {\color{red}0.27476} & 0.14457 & 0.005723 \\ 0 & 0 & {\color{blue}0.19173} & 0.48990 & 0.10794 \\ 0 & 0 & {\color{blue}0.00833} & 0.10833 & 0.05\\ \end{bmatrix} \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 5.17983 \\ 1.17265 \\ 0.95685 \\ 1 \end{bmatrix} \)
Skalujemy teraz wiersz trzeci by dostać jedynkę na przekątnej
\( 3^{rd}=3^{rd}\frac{1}{A_{3,3}}=3^{rd}\frac{1}{0.27476}=3^{rd}3.63953 \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & 1 & 2.99287 & 1.21226 & 0.04676 \\ 0 & 0 & 0.27476{\color{red}*3.63953} & 0.14457{\color{red}*3.63953} & 0.005723{\color{red}*3.63953} \\ 0 & 0 & 0.19173 & 0.48990 & 0.10794 \\ 0 & 0 & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 5.17983 \\ 1.17265{\color{red}*3.63953} \\ 0.95685 \\ 1 \end{bmatrix} \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & 1 & 2.99287 & 1.21226 & 0.04676 \\ 0 & 0 & {\color{red}1} & {\color{red}0.52617} & {\color{red}0.02082} \\ 0 & 0 & 0.19173 & 0.48990 & 0.10794 \\ 0 & 0 & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 5.17999 \\ {\color{red}4.26763} \\ 0.95683 \\ 1 \end{bmatrix} \)
Użyjemy teraz trzeciego wiersza żeby wygenerować zera w trzeciej kolumnie pod przekątną.
\( 4^{th}=4^{th}-3^{rd}A_{4,3}=4^{th}-3^{rd}0.19173 \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & 1 & 2.99287 & 1.21226 & 0.04676 \\ 0 & 0 & 1 & 0.52617 & 0.02082 \\ 0 & 0 & {\color{red}0} & 0.48990{\color{red}-0.52618*0.19173} & 0.10794{\color{red}-0.02083*0.19173} \\ 0 & 0 & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \\ \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 5.17983 \\ 4.26789 \\ 0.95685{\color{red}-4.26789*0.19173} \\ 1 \end{bmatrix} \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & 1 & 2.99287 & 1.21226 & 0.04676 \\ 0 & 0 & 1 & 0.52617 & 0.02082 \\ 0 & 0 & {\color{red}0} & {\color{red}0.38901} & {\color{red}0.10395} \\ 0 & 0 & 0.00833 & 0.10833 & 0.05\\ \end{bmatrix} \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 5.17983 \\ 4.26789 \\ {\color{red}0.13860} \\ 1 \end{bmatrix} \)
\( 5^{th}=5^{th}-3^{rd}A_{5,3}=5^{th}-3^{rd}0.00833 \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & 1 & 2.99287 & 1.21226 & 0.04676 \\ 0 & 0 & 1 & 0.52617 & 0.02082 \\ 0 & 0 & 0 & 0.38901 & 0.10395 \\ 0 & 0 & {\color{red}0} & 0.10833{\color{red}-0.52618*0.00833} & 0.05{\color{red}-0.02082*0.00833}\\ \end{bmatrix} \\ \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 5.17983 \\ 4.26789 \\ 0.13860 \\ 1{\color{red}-4.26789*0.00833} \end{bmatrix} \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & 1 & 2.99287 & 1.21226 & 0.04676 \\ 0 & 0 & 1 & 0.52617 & 0.02082 \\ 0 & 0 & 0 & 0.38901 & 0.10395 \\ 0 & 0 & 0 & {\color{red}0.10395} & {\color{red}0.04983}\\ \end{bmatrix} \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 5.17983 \\ 4.26789 \\ 0.13860 \\ {\color{red}0.96445} \end{bmatrix} \)
Następnie przeskalujemy czwarty wiersz tak żeby miał jedynkę na przekątnej
\( 4^{th}=4^{th}\frac{1}{A_{4,4}}=4^{th}\frac{1}{0.38901}=4^{th}2.5706 \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & 1 & 2.99287 & 1.21226 & 0.04676 \\ 0 & 0 & 1 & 0.52617 & 0.02082 \\ 0 & 0 & 0 & 0.38365{\color{red}*2.5706} & 0.10395{\color{red}*2.5706} \\ 0 & 0 & 0 & 0.10395 & 0.04982\\ \end{bmatrix} \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 5.17983 \\ 4.26789 \\ {\color{red}0.13860*2.5706} \\ 0.91110 \end{bmatrix} \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & 1 & 2.99287 & 1.21226 & 0.04676 \\ 0 & 0 & 1 & 0.52617 & 0.02082 \\ 0 & 0 & 0 & {\color{red}1} & {\color{red}0.26721} \\ 0 & 0 & 0 & 0.10395 & 0.04983 \\ \end{bmatrix} \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 5.17983 \\ 4.26789 \\ {\color{red}0.35629} \\ 0.96444 \end{bmatrix} \)
Użyjemy teraz czwarty wiersz żeby wygenerować zero w czwartej kolumnie pod przekątną
\( 5^{th}=5^{th}-4^{rd}A_{5,4}=5^{th}-4^{rd}0.10395 \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & 1 & 2.99287 & 1.21226 & 0.04676 \\ 0 & 0 & 1 & 0.52617 & 0.02082 \\ 0 & 0 & 0 & 1 & 0.26721 \\ 0 & 0 & 0 & {\color{red}0} & 0.04983{\color{red}-0.26721*0.10395} \\ \end{bmatrix} \\ \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 5.17983 \\ 4.26789 \\ 0.35629\\ 0.96444{\color{red}-0.35629*0.10395} \end{bmatrix} \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & 1 & 2.99287 & 1.21226 & 0.04676 \\ 0 & 0 & 1 & 0.52617 & 0.02082 \\ 0 & 0 & 0 & 1 & 0.26721 \\ 0 & 0 & 0 & {\color{red}0} & {\color{red}0.022053} \\ \end{bmatrix} \nonumber\\ \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 5.17983 \\ 4.26789 \\ 0.35629 \\ {\color{red}0.92740} \end{bmatrix} \)
Na końcu przeskalujemy piąty wiersz tak żeby miał jedynkę na przekątnej
\( 5^{th}=5^{th}\frac{1}{A_{5,5}}=5^{th}\frac{1}{0.022053}=5^{th}45.351 \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & 1 & 2.99287 & 1.21226 & 0.04676 \\ 0 & 0 & 1 & 0.52617 & 0.02082 \\ 0 & 0 & 0 & 1 & 0.26721 \\ 0 & 0 & 0 & {\color{red}0} & {\color{red}1} \\ \end{bmatrix} \\ \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 5.17983 \\ 4.26789 \\ 0.35629\\ 0.92740{\color{red}*45.351} \end{bmatrix} \)
\( \begin{bmatrix} 1 & 4.61553 & 2 & 0.07693 & 0 \\ 0 & 1 & 2.99287 & 1.21226 & 0.04676 \\ 0 & 0 & 1 & 0.52617 & 0.02082 \\ 0 & 0 & 0 & 1 & 0.26721 \\ 0 & 0 & 0 & {\color{red}0} & {\color{red}1} \\ \end{bmatrix} \\ \begin{bmatrix} u_{1} \\ u_{2} \\ u_{3} \\ u_{4} \\ u_{5} \\ \end{bmatrix} = \begin{bmatrix} 9.23105 \\ 5.17983 \\ 4.26789 \\ 0.35629\\ {\color{red}42.059} \end{bmatrix} \)
Rozpisujemy \( 1u_1+4.61553u_2+2u_3+0.07693u_4=9.23105 \\ 1u_2+2.99287u_3+1.21226u_4+0.04676u_5=5.17983 \\ 1u_3+0.52617u_4+0.02082u_5=4.26789 \\ 1u_4+0.26721u_5=0.35629 \\ 1u_5=42.059 \)
I rozwiązujemy od ostatniego równania
\( u_5={\color{blue}42.059} \\ u_4=0.35629-0.26721u_5 =0.35629-0.26721*{\color{blue}42.059}={\color{red}-10.88229} \\ u_3=4.26789-0.52617u_4-0.02082u_5 =4.26789-0.52617{\color{red}-10.88229}-0.02082{\color{blue}42.059}={\color{green}9.11815} \\ u_2=5.17983-2.99287u_3-1.21226u_4-0.04676u_5= \\ =5.17983-2.99287{\color{green}9.11815}-1.21226{\color{red}-10.88229}-0.04676{\color{blue}42.059} = {\color{cyan}-10,88412} \\ u_1=9.23105-4.61553u_2-2u_3-0.07693u_4 =\\=9.23105-4.61553{\color{cyan}-10,88412}-2{\color{green}9.11815}-0.07693{\color{red}-10.88229} = 42.0679 \)
Tym razem nasze rozwiązanie jest poprawne.

Algorytm elimiancji Gaussa należy wzbogacić o zamianę wierszy


1 for irow=1,N //pętla po wierszach
2 Znajdź największą wartość w kolumnie A(row:N,irow)
3 Oznacz wiersz z największą wartością irowmax
4 Zamień wiersze irow, irowmax
5 for icol=irow+1,N // skalowanie wiersza
6 A(irow,icol)=A(irow,icol)/A(irow,irow)
7 end loop
8 b(irow)=b(irow)/A(irow,irow)
9 A(irow,irow)=1.0
10 for irow2=irow+1,N //odejmowanie wierszy
11 for icol=irow2+1,N
12 A(irow2,icol) = A(irow2,icol) - A(irow2,irow) * A(irow,icol)
13 end loop
14 b(irow2)=b(irow2)-A(irow2,irow)*b(irow2)
15 A(irow2,irow)=0
16 end loop
17end loop


Ostatnio zmieniona Sobota 09 z Lipiec, 2022 08:19:36 UTC Autor: Maciej Paszynski
Zaloguj się/Zarejestruj w OPEN AGH e-podręczniki
Czy masz już hasło?

Hasło powinno mieć przynajmniej 8 znaków, litery i cyfry oraz co najmniej jeden znak specjalny.

Przypominanie hasła

Wprowadź swój adres e-mail, abyśmy mogli przesłać Ci informację o nowym haśle.
Dziękujemy za rejestrację!
Na wskazany w rejestracji adres został wysłany e-mail z linkiem aktywacyjnym.
Wprowadzone hasło/login są błędne.